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Abstract 

For most physical processes the tracks observed in the muon sta- 
tions must be matched with the corresponding tracks in the inner 
tracker, the external muon system providing muon identification and 
triggering but the tracker points giving the precise momentum mea- 
surement at lower momenta. For high momenta the momentum res- 
olution is much improved by the interconnection of inner and outer 
measurements. The matching of outer and inner measurements is more 
delicate in case of muons embedded in jets. A study of the matching 
procedure was carried out using samples of bb jets at high Pt, requiring 
bb ^ fi decays. 



1 Motivation 



Identification and measurement of muons in the inner tracker using its higli 
spatial resolution and excellent momentum measurement in the 4T field is a 
very important and demanding task, particularly for muons in jets. In the 
muon system the occupancy is low and pattern recognition is relatively easy, 
provided the particle reaches the muon system at all, which is the case for 
the barrel is > 4 GcV. For non-isolated muons the main difficulty is to 
match the track observed in the muon stations with the correct candidate 
in the inner tracker in case of the high track density inside a jet. These 
situations arise for example in muon tagging of 6-jets, in 6 — > J/il){-^ jx 
^) + X cascades in h and top physics, in SUSY(^, g) physics and most 
dramatically in very high (TeV scale) jets such as heavy flavour (5, c) 
QCD jets giving rise to "muon bundles" within jets. 

The matching strategy is based on the use of the full parameter space, 
i.e. on both the reconstructed momentum and the position in a given plane 
for the two tracks, one in the muon system and the other in the inner 
tracker. The quantities monitored are the efficiency of matching (number of 
muons matched over number of muons reconstructed), the muon matching 
impurity (1 - number of muons matched over number of tracks matched) 
and generation of "ghost" and "fake" tracks. 

2 CMS model and data sample 

The results discussed here are obtained using a detailed GEANT description 
of the complete CMS detector of TDR design. In all simulations the effects 
of the 4T field are fully included. 

We have performed a GEANT simulation of the response to bb jets with 
transverse momenta P^*^* of 50, 100 and 200 GeV (which corresponds to the 
Pj range of 6-jets in ti production) and at least one muon inside each jet 
with P^ > 5 GeV in the rapidity region -1.1 < 77 < 1.1. 1000 events were 
generated and reconstructed in the barrel part at each energy point. 

The results are obtained using the standard CMS reconstruction proce- 
dure (data card 'MUON' 2 and 'TRAK' 3) existing in the CMSIM version 
115. 

Pattern recognition in the muon stations is done by building track seg- 
ments from the reconstructed hits separately in the (r,(/)) and (r,z) planes; 
(r^,rz) track segments in the same muon station and in the same or contigu- 
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ous (r,0) sectors are considered to form a pair used to perform the Kalman 
filter. This is done starting from the outermost stations and proceeding 
towards the inner stations with a constrained fit to the impact point (IP) 
with default values of cr^;^ = 15 fim and = 5.3 cm. Since the decay 
mean transverse length (impact parameter) of the 6-quarks is about 500 nm 
another reconstruction was made using a constrained fit to the impact point 
with = 500 nm and a^^ = 5.3 cm. 

The Global Track Finder (GTF) was used for pattern recognition in 
the CMS inner tracker. The GTF algorithm adopted for the CMS tracker 
must process a large number of hits per event: typically 5.10^ hits at low 
luminosity and ten times more at high luminosity. To overcome the severe 
combinatorial problems, the concept of road preselection is used in the first 
stage of the algorithm. In the second stage a Kalman filter is used to carry 
out final hit selection and track fitting. 

3 Track matching 

Given two track segments, one in the muon system and the other in the in- 
ner tracker, one can extrapolate both segments to a common plane or point, 
check the difference in momentum and position and select the best combi- 
nation among the candidates. The correct combination is known from the 
simulation procedure and it is therefore possible to calculate the efficiency 
of matching. We use as a matching criterion the minimal distance computed 
from 77, and/or Pp*^ of the two track candidates at the outer plane (TRCI 
- CMSIM convention) of the inner tracker and at the impact point. As- 
sociation of reconstructed tracks to GEANT ones is done on the hit basis. 
The reconstructed track is associated to a given GEANT particle number 
if more than half of the hits assigned to the track are from that particular 
GEANT track. If there is no such dominating track, the reconstructed track 
is called a "ghost". If more than one reconstructed track is matched to a 
given GEANT number it is called a " fake" . 

The 77 dependence of the reconstruction efficiency obtained for the barrel 
muon station, shown on fig. la) is rather uniform except for the region 
between the central and next-to-central rings (| r/ | < 0.2) and for | rj \ 
> 0.9. To avoid the infiuence of the r] dependence of the reconstruction 
efficiency on the matching results the tracks with | r/ | > 0.9 are removed. 
The P( reconstruction efficiency (fig. lb)) is more than 95% for the tracks 
with P^ > 10 GeV. In fig. Ic) is given the "type" of reconstructed tracks in 



2 



the muon stations (solid line) for which must be found the best candidate 
between the tracks reconstructed in the tracker (error bars). The comparison 
between reconstructed and generated variables: Pt, r] and (f) is given in fig. 
Id) - fig. li) both for the barrel muon stations and for the inner tracker, for 
P^^* = 100 GeV. We find a Pj muon resolution of about 10% for the muon 
stations alone, and of < 1% for the inner tracker, in good agreement with 
the values given in the Muon and Tracker TDR's/^'^/. 

To estimate the matching quality a is formed from the differences 
between reconstructed r], (p and Pj in the two detectors for all reconstructed 
track combinations at the outer tracker plane and at the impact point: 



Including in the x^, equation (1), the contribution from PJ'^'^ in the muon 
stations and the tracker improves the selection of the matching tracks mainly 
for the high transverse momenta. A minimal of the matching tracks has 
been used as a matching criterion. 

The Pt dependence of the matching efficiency is shown on fig. 2 for P^^* 
= 50, 100 and 200 GeV at the outer plane TRCI of the inner tracker and 
at the interaction point with the two assumptions for the constrained fit to 
the impact point with ax,y = 15 /xm and ax,y = 500 fim. The degradation 
of the matching efficiency at higher P^^* is is due to the increase of charged 
particle multiplicity inside jets which makes more difficult the correct asso- 
ciation of the tracks observed in the muon stations with the track candidates 
in the inner tracker. The matching efficiency shows less variation when a 
constrained fit at the impact point is performed using ax,y = 500 /xm. 

Results obtained from the matching of muons inside the bb jets at P|*^* 
of 50, 100 and 200 GeV are summarized in Table 1. The "type" of the 
matched tracks indicates a higher purity at the impact point and at lower 
jet momenta. For the higher multiplicity and more collimated jets at P^^* 
= 200 GeV, the impurity in matching degrades to a ~ 7% level. Rejecting 
tracks with an increasing momentum threshold < P^"*, a high matching 
efficiency and purity of track matching can be obtained for high Pj^*, but at 
the expense of a substantial reduction in statistics. The number of "ghost" 
tracks generated at reconstruction in the muon stations is negligible. The 
number of "fake" tracks shows some tendency to increase with jets momenta, 
and is at level of 1% - 2%. 
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After performing the matching procedure we can assign to the tracks 
observed in the muon station the precise measurements from the tracker. 
The distributions of the differences between reconstructed and generated 
vahies of P^, 77 and (p at the impact point for the matched tracks at P^^* = 
100 GcV are given in fig. 3. These distributions are practically the same as 
in fig. 2, where we compared reconstructed and generated variables for the 
inner tracker, confirming the good quality of the matching criterion. The 
bottom row in fig. 3 gives the "type" of the tracks reconstructed in the 
muon stations and the "type" of matched tracks at the impact point. 

The quality of track matching can be measured by the value of the "pull" : 

p ^trPtr ^msPms |'2\ 

where: 

zjj., 'Z'ms is the track charge in the tracker and the muon station; 

Ptr, Pms is the reconstructed momentum in tracker and muon station; 

cov(l/p.l/p) is the covariancc matrix. 

A Gaussian fit to the pull shown at the bottom of fig. 3 gives cr ~ 1.0 
indicating a good accuracy for the selection and matching of the tracks. 



4 Summary and conclusion 

The main results concerning the matching of muon tracks in (65) jets between 
the muon stations and inner tracker are: 

Using a simple minimal distance criterion between ry, ^ and Pj of in- 
dependent track measurements at the impact point one can obtain a high 
matching efficiency with a contamination of not more than 4% for Pf * of 
50 GeV and 6% for 100 GeV. For the more coUimated jets with Pf * = 200 
GeV, the impurity in matching is at a ~ 7% level. 

The matching efficiency shows less variation when a constrained fit at 
the impact point is performed with ax,y = 500 jim. To further improve 
the matching efficiency and impurity a more sophisticated software than 
currently available must be used, for example a constrain fir to the secondary 
hh jets vertices obtained from the inner tracker can be performed. 

Rejecting tracks with an increasing momentum cut-off < Pf"*, a high 
purity of track matching can be obtained for hard jets, but at the expense 
of a substantial reduction of the statistics. 
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Table 1: Track matching results. 



P^^' = 50 GeV 


TRCI 


IP <^x,y = 15 


IP <^x,y = 500 /xm 


PJ"' (GcV) 


loss of stat. (%) 


efficiency (%) 


impurity (%) 


efficiency (%) 


impurity (%) 


efficiency (%) 


impurity (%) 
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96.42 


3.64 ± 0.47 


96.22 


4.08 ± 0.50 


96.83 


3.23 ± 0.45 


6 


2.68 


96.85 


2.89 ± 0.43 


96.98 


3.02 ± 0.48 


97.06 


2.81 ± 0.43 


7 


7.07 


97.53 


2.00 ± 0.37 


97.73 


1.86 ± 0.33 


97.63 


1.89 ± 0.36 


8 


11.98 


97.97 


1.60 ± 0.34 


98.26 


1.38 ± 0.30 


98.16 


1.11 ± 0.29 


9 


17.97 


98.13 


1.41 ± 0.33 


98.52 


1.09 ± 0.24 


98.43 


0.79 ± 0.25 


10 


24.79 


98.47 


0.85 ± 0.27 


98.73 


0.68 ± 0.23 


98.46 


0.60 ± 0.23 


"ghost" tracks(%) 


0.13 = 


b 0.09 


0.06 ± 0.06 


"fake" tracks(%) 


1.47 ± 0.30 


1.45 ± 0.39 


Pj«^* = 100 GeV 


TRCI 


IP <7a;,y = 15 /im 


IP <^x,y = 500 


P™' (GcV) 


loss of Stat. (%) 


efficiency (%) 


impurity (%) 


efficiency (%) 


impurity {%) 


efficiency (%) 


impurity (%) 







93.57 


6.55 ± 0.61 


94.40 


6.24 ± 0.62 


94.87 


5.89 ± 0.59 


6 


2.22 


93.93 


6.13 ± 0.60 


95.18 


5.24 ± 0.56 


95.61 


5.05 ± 0.55 


7 


4.32 


94.38 


5.68 ± 0.59 


95.66 


4.71 ± 0.58 


96.43 


4.13 ± 0.51 


8 


6.49 


94.78 


5.28 ± 0.57 


96.16 


4.16 ± 0.54 


96.94 


3.44 ± 0.47 


9 


9.57 


95.15 


4.85 ± 0.56 


96.17 


3.96 ± 0.56 


97.12 


3.15 ± 0.46 


10 


13.03 


95.66 


4.27 ± 0.54 


96.73 


3.27 ± 0.48 


97.30 


2.84 ± 0.44 


15 


26.03 


97.39 


2.36 ± 0.44 


98.06 


1.77 ± 0.33 


98.83 


1.17 ± 0.31 


"ghost" tracks(%) 


0.08 = 


b 0.07 


0.10 ± 0.08 


"fake" tracks(%) 


1.78 ± 0.41 


1.36 ± 0.38 


P^"^ = 200 GeV 


TRCI 


IP '^x,y = 15 


IP '^x,y = 500 pim 


PJ"' (GcV) 


loss of stat. (%) 


efficiency (%) 


impurity (%) 


efficiency (%) 


impurity (%) 


efficiency (%) 


impurity (%) 


5 





90.17 


9.99 ± 0.74 


92.55 


8.07 ± 0.67 


94.12 


6.90 ± 0.62 


6 


0.60 


90.60 


9.51 ± 0.72 


93.00 


7.51 ± 0.69 


94.33 


6.59 ± 0.61 


7 


1.62 


91.18 


8.94 ± 0.71 


93.42 


6.98 ± 0.66 


94.70 


6.12 ± 0.60 


8 


3.25 


91.72 


8.40 ± 0.69 


93.69 


6.72 ± 0.67 


95.00 


5.48 ± 0.57 


9 


5.05 


91.94 


8.17 ± 0.69 


94.15 


6.14 ± 0.68 


95.24 


5.07 ± 0.56 


10 


6.80 


92.24 


7.88 ± 0.68 


94.56 


5.75 ± 0.52 


95.21 


5.04 ± 0.56 


15 


17.63 


93.69 


6.38 ± 0.65 


95.98 


4.23 ± 0.53 


96.44 


3.70 ± 0.51 


20 


23.95 


94.86 


5.07 ± 0.62 


96.75 


3.25 ± 0.41 


97.50 


2.65 ± 0.46 


"ghost 


' trackB(%) 


0.18 : 


h 0.10 


0.18 ± 0.10 


"fake" tracks(%) 


1.99 ± 0.34 


2.11 ± 0.35 
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Figure 1: r) (a) and (b) dependence of the reconstructed efficiency, "type" 
of reconstructed tracks (c) in the muon stations (sohd hne) and in the tracker 
(error bars) and comparison between reconstructed and generated variables 
Pt, 7] and (p in the barrel muon station (d - f) and in the inner tracker (g - 
i) for Pf * = 100 GeV. 
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Figure 2: P( dependence of the matching efficiency for P^^ = 50, 100 and 
200 GeV at the outer plane TRCI of the inner tracker and at the impact 
point with two assumptions for the constrained fit to the impact point with 
(^x,y = 15 /xm and a^^y = 500 /xm respectively. 
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Figure 3: Comparison of reconstructed and generated values of Pj, rj and cp, 
"type" and "pull" for matched tracks at impact point for P^^* = 100 GeV. 
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